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Abstract 

We consider possible leptonic three-body decays of spin-1/2, charge-asymmetric dark matter. 
Assuming a general Dirac structure for the four-fermion contact interactions of interest, we study 
the cosmic-ray electron and positron spectra and show that good fits to the current data can be 
obtained for both charged-lepton-flavor-conserving and flavor-violating decay channels. We find 
that different choices for the Dirac structure of the underlying decay operator can be significantly 
compensated by different choices for the dark matter mass and lifetime. The decay modes we 
consider provide differing predictions for the cosmic-ray positron fraction at energies higher than 
those currently probed at the PAMELA experiment; these predictions might be tested at cosmic-ray 
detectors like AMS-02. 
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I. INTRODUCTION 



Cosmic rays have been studied extensively at various earth-, balloon- and satellite-based 
experiments. Recently, the PAMELA satellite has observed an unexpected rise in the cosmic- 
ray positron fraction from approximately 7 to 100 GeV This feature is not explained by 
the expected background from the secondary production of cosmic-rays positrons. Moreover, 
observations of the total flux of electrons and positrons by Fermi-LAT ^ and H.E.S.S. jsl 
also show an excess over the predicted background, up to an energy of ~ 1 TeV. The presence 
of nearby pulsars could provide an astrophysical explanation for these observations 
Nevertheless, more exotic scenarios remain possible. The annihilation of dark matter in 
the galactic halo to electrons and positrons provides one such possibility, though generic 
annihilation cross sections must be enhanced by a large boost factor in order to describe the 
data js, 7|. Alternatively, the excess could be explained by a TeV-scale decaying dark matter 
candidate. (For example, see Ref. jsj; for a recent review, see Ref. jol.) In this scenario, fits 
to the cosmic-ray data indicate that the dark matter must decay primarily to leptons with 
a lifetime of O(1026) s. 

While the thermal freeze-out of weakly-interacting, electroweak-scale dark matter can 
naturally lead to the desired relic density, this is not the only possible framework that can 
account for the present dark matter abundance. Recently proposed asymmetric dark matter 
models relate the baryon or lepton number densities to the dark matter number density. 



motivated by the fact that these quantities are not wildly dissimilar 10l-ll3|. 



asymmetric dark matter models have been constructed, for example, in Refs. 
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eV-scale 



13|. The 



asymmetry between dark matter particles and antiparticles can lead to differences in the 
primary cosmic-ray spectra of electrons and positrons, with potentially measurable conse- 



quences Q, |15|. Evidence for such charge asymmetric dark matter decays would disfavor 



n 



the pulsar explanation of the excess 



15|. In addition, charge asymmetric dark matter 



decays may allow one to discern whether dark matter decays are lepton-flavor- violating 16 |. 
For example, the cosmic-ray spectra that one expects if dark matter decays symmetrically 
to e^/i~ and e~fi'^ are indistinguishable from those obtained by assuming flavor-conserving 
decays to e^e~ and fi~^fi~ with equal branching fraction; the same is not true if the dark 
matter decays asymmetrically to e~^fi~ alone, 100% of the time. 



Refs. 



15| and 16| study the cosmic-ray spectra assuming a number of two-body charge- 
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asymmetric dark matter decays, with the latter work focusing on lepton-flavor-violating 
modes. In this paper, we extend this body of work to charge-asymmetric three-body decays 
and, in particular, to modes that violate lepton flavor. We assume a spin-1/2 dark matter 
candidate that decays via four-fermion contact interactions to two charged leptons and a 
light, stable neutral particle. For the present purposes, the latter could either be a standard 
model neutrino or a lighter dark matter component. Four-fermion interactions have a long 
history in the development of the weak interactions, and one can easily imagine that dark 
matter decays could be the consequence of operators of this form, generated by higher-scale 
physics. Moreover, the possible presence of a neutrino in the primary decay may lead to 
interesting signals at neutrino telescopes [l^]. Unlike the two-body decays already considered 
in the literature, the precise energy distribution of the decay products is affected by the 
Dirac matrix structure of these contact interactions, which is not known (unless a model is 
specified). By considering the most general possibilities, we show that different choices for 
the Dirac structure of the decay operators defined in Sec. [TTlcan be substantially compensated 
by different choices for the dark matter mass and lifetime r^; while the best fit values of 
these parameters change, the predicted spectra are not dramatically altered. On the other 
hand, we find that the flavor structure of the decay operator has a more significant effect. 
Assuming various lepton-flavor-conserving and flavor-violating decay modes, we compute 
the resulting cosmic-ray spectra, performing fits to the data to determine the optimal 
dark matter masses and lifetimes. Like Refs. [isl, [l^ , we obtain predictions for these spectra 
at energies that are higher than those than can be probed accurately now. Future data 
from experiments like AMS-02 may provide the opportunity to test these predictions, 
and evaluate them relative to other interpretations of the cosmic-ray positron excess. 

This letter is organized as follows. In the next section, we discuss the assumed form of 
the dark matter operators. In Sec. Illlt we present the results of our numerical analysis and 
in Sec. \IV\ we discuss our results and directions for future work. 



II. FOUR-FERMION OPERATORS 

We consider a spin-1/2 dark matter candidate tp that decays to d-fd-Jv where i and j are 
generation indices and v represents a light, neutral particle. We assume that v is either 
a standard model neutrino or a secondary dark matter component that is much lighter 
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than ip and contributes negligibly to the relic density. In the present analysis, the exact 
nature of the light neutral state will be irrelevant since its effect on our results will come 
solely from kinematics. We focus on the simplest scenario, in which there are no additional 
decay channels involving the charge conjugate of z/, and consider the possible four-fermion 
operators that contribute to the decays of interest. We work directly with the operators that 
may appear after the standard model electroweak gauge symmetry is spontaneously broken; 
for any operator found to have phenomenologically desirable properties, one may easily 
construct a gauge-invariant origin after the fact. Note that the production of a neutrino in 
the primary decay may have interesting phenomenological consequences (see, for example. 



Ref. [l7j), which provides a separate motivation for our three-fermion final state. Once this 
choice is made, the dark matter spin must be 1/2 if the underlying theory is renormalizable ^. 

The problem of parametrizing an unknown decay amplitude of one spin-1/2 particle to 
three distinct spin-1/2 decay products was encountered in the study of muon decay, before 
the standard model was well established. The most general decay amplitude M. can be 



parametrized by |20| 

iM = igY^ [u{j>o)OiU^] \u{j>^)Oi{c, + c^7^)t;(p+)] , (2.1) 

i 

where p± and po are the momenta of the decay products, labeled according to their electric 
charge, and the Oj, z = 1 ■ ■ ■ 5 are elements of the set of linearly independent matrices 

o = {l,7^a^^7^^7'}• (2.2) 

The Ci and c- are complex coefficients. Terms involving the contraction of spinor indices 
that link different pairs of spinor wave functions can be recast in the form of Eq. (12. ip via 
Fierz transformations. Since the final state particles are much lighter than the dark matter 
candidate (which is at the TeV scale), we can safely neglect their masses. 

Since the neutral final state particle is stable, the energy spectra of electrons and positrons 
that are observed at cosmic-ray observatories are determined by the energy spectra of the 
the charged leptons, ^^ and that are produced in the primary decay; this follows from 
the differential decay distribution 

^ ^ -d-MP), (2.3) 



r dEodE± GAn^m, 



^ For a model with flavor-conserving, three-body decays involving a final-state gravitino, see Ref. [l9| . 



where (|A4p) is the spin-summed/ averaKed squared amphtude. We evaluate this quantity 



exactly from Eq. (12. ip using FeynCalc 2l|, and compute the £^ energy distribution by inte- 
grating over the neutral lepton energy Eq. We find that the result contains terms quadratic 
and cubic in E±; however, since the distribution must be normalized to unity, the result has 
the following simple parametrization: 



1 dT IE. 



r dE. - -2 



e± + 64 - -e± . 

6 J 



E. 



± 



(2.4) 



The requirement that this expression remains positive over the kinematically accessible range 
< E± < restricts the parameters {_|_ and to fall within the range 

< C± < 96 . (2.5) 

The ^± are generally complicated functions of the operator coefficients q and we provide 
these in the appendix. In the present analysis, however, the exact relations are not par- 
ticularly important; by leaving and as fitting parameters, one obtains very similar 
predicted spectra, independent of the choice of the The fact that some solution exists 
for any desired Dirac structure of the underlying four-fermion operator makes it potentially 
easier to construct explicit models. Though we reserve the task of model-building to future 
work, it is worth noting, for example, that the operator 

qRR ^ + + ^^y^ , (2.6) 

corresponding to = 96 and = 48, is a particularly interesting choice, since it is already 
gauge invariant under the standard model gauge group and may provide a simple starting 
point for constructing a plausible ultraviolet completion. 



We computed the electron and positron spectra using PYTHIA [22|, taking into account 
the energy distributions of the primary leptons and l~ . As a cross check, we have written 
code that incorporates Eq. (12. 3p . computed directly from a choice of the underlying four- 
fermion operator, as well as code that incorporates only the distributions Eq. (12. 4p . for the 
corresponding values of ^+ and We have also compared output from different versions 
of our code, based on PYTHIA 6.4 and PYTHIA 8.1, respectively^. Results from these 
different approaches were found to be agreement. 



^ Note that PYTHIA 6.4 does not automatically take into aceount neutron decay, which we include by 
modifying the program's decay table. 
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III. COSMIC-RAY SPECTRA 



To compute the relevant cosmic-ray fluxes, one must take into account that electrons and 
positrons produced in dark matter decays must propagate through the galaxy before reaching 
earth. While modeling this propagation is now standard in the literature on decaying dark 
matter scenarios, we briefly summarize our approach so that our discussion is self contained 
and our assumptions are manifest. 



A. Cosmic-Ray Propagation 

Let r be a position with respect to the center of the Milky Way Galaxy. We assume the 



spherically symmetric Navarro-Frenk- White dark matter halo density profile 



23| 



^3 

p{r) = Po I ,^ X2 ^ (3.1) 

where po — 0.26 GeV/cm^ and ~ 20 kpc. The production rate of electrons/positrons per 
unit energy and per unit volume is then given by 

Q(E,r)^^-^(L^M\, (3,2) 



dE 

where and are the dark matter mass and lifetime, respectively, and dNe±/dE is the 
energy spectrum of electrons/positrons produced in the dark matter decay. Let fe±{E,r) 
be the number density of electrons/positrons per unit energy. Then, fe±{E,r) satisfies the 



transport equation 



24| 



= i^(i?)VVe±(^,r) + A [b{E)U{E,r)] + QiE,r). (3.3) 



We assume the MED propagation model described in Ref. |25|] for which 

K{E) = 0.0112e°-™kpcVMyr (3.4) 

and 

6(E) = lO^VGeV/s, (3.5) 

where e = -E/(lGeV). The diffusion zone is approximated as a cylinder with half- height 
L = 4 kpc and radius i? = 20 kpc. We require fe±{E, r) to vanish at the boundary of this 
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zone. The solution at the hehospheric boundary is then given by 26 1 



U{E) = ^ / dE'G^.{E,E')'^J^^. (3.6) 
m^r^ J dE' 







The Green's function, Ge±{E, E'), can be found in Ref. 26j. The interstellar flux of elec- 



trons/positrons created in dark matter decays is then given by 

^e^iE) = ^fe4E), (3.7) 

where c is the speed of light. 

For the background fluxes, we assume the Model proposed by the Fermi collabora- 
tion |2j 



28 



and 



= ( 1 + 0224^'^ ' GeV-'m-s-'=r-' (3.8) 
**(^' = f 1 + 0,0002'' « + 2^-*'"' ' G^V-'m-s-'sr- , (3.9) 



where, as before, e = E/{1 GeV). 

At the top of the Earth's atmosphere, these fluxes must be corrected to account for the 



effects of solar modulation [28 
the interstellar (IS) flux by 



. The flux at the top of the atmosphere (TOA) is related to 



C=?^(i?TOA) = ^e(i^is), (3.10) 

^IS 



where Eis = -Etoa + and \e\4>F = 550 MeV. 
The total electron-positron flux is given by 

^tot _ $DM^^^ ^ $°+^(^) + k^l^.^{E) + ^1+^{E) , (3.11) 

where k is a free parameter which determines the normalization of the background electron 
flux. In our numerical analysis, we find that the best fit values of k never deviate by 
more that two percent from 0.84 and that fixing k at this value has a negligible effect on 
the goodness of fits and our predicted spectra. Therefore, we set k = 0.84 henceforth to 
reproduce the cosmic-ray spectra at low energies. The positron fraction is given by 

PP(.) . «t^. ,3,:2) 

e 
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B. Results 



In the propagation model described above, the only remaining undetermined quantities 
are m^, r^, dN^+jdE and dNf,-/dE. The electron and positron energy spectra, dN^+fdE 
and dN^-fdE, are determined by and by a set of parameters which we describe in the 
following paragraph. 

We consider dark matter decays of the form — )■ (f^J^ where if is a charged lepton of 
the i^^ generation. There are nine such decay channels, and we require 

J2B{iiiJu) = l, (3.13) 

id 

where the B{ifi~v) are branching fractions. For decays involving more than one channel, 

. ,3.4) 

where {dNe±/dE)-j is the electron/positron energy spectrum for ip — )■ ifiji^- In Sec. [Tll 
we showed that the energy spectra of the charged leptons in the decay — )■ if^J^ 
characterized by the ordered pair (,^+,,^„), where < < 96. We also showed that 
{dNe±/dE)-j is entirely determined by and (^+,^_). For decays involving more than one 
decay channel (e.g., ip — e^^~v and ip — fi'^r^u), we assume a constant (,^4.,,^_.). Then, 
since the branching fractions are subject to Eq. f l3.13p . we can determine dN^+/dE and 
dNg-/dE by specifying m^, and eight of the nine branching fractions. 

To summarize, when we use the cosmic-ray propagation model described in the previous 
subsection, the resulting positron fraction and total electron-positron flux measured at the 
top of the Earth's atmosphere are determined by 12 parameters: m^, r^, ^_ and eight 
of the nine branching fractions. 

For each of the decay scenarios considered below, we fixed (C+j'C-) ^-nd the branching 
fractions and then performed a fit to the PAMELA, Fermi LAT, H.E.S.S. 2008 and 
H.E.S.S. 2009 data with and as fitting parameters. We allowed to vary in incre- 
ments of 500 GeV, and we allowed r^, to vary in increments of 0.1 x 10^^ s. We consider 
the range E' > 10 GeV, where the effects of a TeV-scale dark matter candidate are rele- 
vant. Where the high-energy and low-energy Fermi data overlap, we have plotted only the 
high-energy data. (We omit from our figures the H.E.S.S. bands of systematic uncertainty.) 
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FIG. 1: The envelope of possible cosmic-ray spectra for ijj — )• t~^t v. Ranges of the fit parameters 
are given in the text. 

Leaving and as free variables, we find that our results are relatively insensitive 
to the choice of (^+,^_). This is demonstrated for the pure decay — t- t^t~v in Fig. [T] 
where we show the envelope of possible cosmic-ray spectra; that is, when we sample the 
(,^+,,^_) parameter space, we find that all of the resulting curves fall between those plotted 
in Fig.[Tl For the example shown, varies between 6.5 and 8.5 TeV while varies between 
0.5 X 10^^ s and 0.7 x 10^^ s; the P^r degree of freedom (x^/d.o.f.) remains between 0.5 
and 0.6. We performed the same analysis on the other decay scenarios discussed below and 
found a similar behavior. As such, we take (^+,^_) = (48,48) for the remaining results that 
we present. 

As a starting point, we show the cosmic-ray spectra for some charged- lepton- flavor- 
conserving decays in Fig. [2l We consider the pure decays il) — ?■ fi^fi^u and — ?■ r^r^u, 
and we also consider the flavor-democratic decay for which B^ifi^u) = 1/3 for all i. 
For ip — fi~^fi~h', we have a x^/d.o.f. of approximately 0.9. For ip t^t^u, we have 
X^/d.o.f. ~ 0.6. And for the flavor-democratic — we have x^/d.o.f. ~ 0.8. These 
are to be contrasted with the flavor- violating decays of Fig. [3l 

We consider three classes of flavor-violating decays: 

^jJ e^/i^z/, ijj — )■ e^T^u, and ijj fi^r^u. (3.15) 

Each class contains two decay channels (e.g., ^p — )■ e"'"/i~z/ and ip e~fi~^i'). We consider all 
six of the pure decays, i.e., decays involving only one channel. We also consider mixtures of 
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FIG. 2: Positron fraction and total electron-positron flux for some charged-lepton-flavor-conserving 
decays. Best fits are shown, corresponding to the following masses and lifetimes: for ^ 

= 3.5 TeV and = 1.5 x lO^^ s; for T+T'iy, = 7.5 TeV and r^, = 0.6 x 10^^ s; for 
the flavor-democratic decay ^jJ t^l~v, = 2.5 TeV and = 1.9 x 10^^ s. 

decay channels belonging to the same class; some representative choices are shown in Fig. |3l 
Note that, for fixed m^, and r^, the total electron-positron flux - which does not distinguish 
between the two electric charges - is the same for any two decays belonging to the same 
class. For this reason, we require only one plot of the total flux in Fig. [3l We find that the 
is relatively fiat as a function of the branching fraction within each class of decays: over 
the range of possible branching fractions, we find that the x^/d.o.f. varies by no more than 
±10% from 1.2, 1.1 and 0.6, for ip — )■ e^fi^u, ip — e^r^u, and ip — fi^r^u, respectively. 
Different choices for the branching fraction within a given class describe the existing data 
well, but provide different predicted spectra that interpolate between the curves shown. 
Note that the distinctive dip in the ^'^e~v and T^e~v positron fractions around 1 TeV is 
due to the hard electron produced in the initial decay; this greatly enhances the electron to 
positron ratio in the high energy bins, leading to a suppression in the positron fraction for 
fixed total flux. 



IV. DISCUSSION 



The results presented in the previous section show that a variety of possible lepton- 
flavor- violating decay modes for a spin-1/2, charge asymmetric dark matter candidate can 
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FIG. 3: Positron fraction and total electron-positron flux for some charged-lepton-flavor-violating 
decays with various sets of branching fractions. Best fits are shown, corresponding to the following 
masses and lifetimes: for ■0 — )> e^fi^u, = 2.0 TeV and = 2.9 x 10^^ s; for t/j — )• e^r^u, 
= 2.0 TeV and = 2.4 x lO^^ s; for i; /j^t^u, = 4.5 TeV and = 1.0 x lO^^ s. 

describe existing data well, as quantified by the per degree of freedom for the best fits to 
the data. Significantly, the results for the predicted positron fraction differ substantially for 
energies above ~ 100 GeV, the maximum for which the PAMELA experiment is sensitive. 
In some case, more precise measurement of the total electron-postron flux around 1 TeV may 
also provide a means of distinguishing these scenarios. Future data from experiments like 
AMS-02 [l^, which can probe these energy ranges of the predicted spectra, may determine 
whether the possibilities discussed in this letter present viable descriptions of the cosmic-ray 
spectrum. 

In the meantime, the present work suggests a number of directions for further study: In 
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the case where the stable, neutral particle in the final state is a standard model neutrino, 
one could study whether the decays of asymmetric dark matter that we have considered 
could be probed at neutrino observatories like IceCube • One could also study additional 
astrophysical bounds on the scenarios described, for example, from the extragalactic gamma 
ray fiux 15|]. One can also attempt to find preferred forms of the underlying four-fermion 
operators (whose effects were parametrized in the present analysis by by studying the 
simplest and best-motivated models that provide for their origin. Work in these directions 
is in progress and will be described in a longer publication. 



Acknowledgments 

This work was supported by the NSF under Grants PHY-0757481 and PHY-1068008. 
A.C.'s undergraduate research was supported, in part, by a William & Mary Honors Fellow- 
ship. CD.C gratefully acknowledges support from a William & Mary Plumeri Fellowship. 

Appendix A: The Parameters ^± 



The parameters ^± may be expressed in terms of the operator coefficients Cj and c- defined 
in Eq. dSD, 



48 



c'^N± c + c'^N± c' 



(Al) 



cWc + c'^Dd ' 

where c = [ci, C2, C3, C4, c^]'^ and c' = [c[, c'2, C3, C4, c'^]^. The five-by-five matrices N± and D 
are given by 



( I ^2 \ 
6 ±2 
^2 40 ^2 
±2 6 
\ ^2 1 / 



and D 



/l o\ 

4 
24 
4 
\0 1 / 



(A2) 
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